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FIGURES 


Electron  Fluxes  from  the  S3-3  108  through  65^  keV 
Channels  Obtained  on  28  July  1976  While  the  Satellite 
Was  below  the  Stable  Trapping  Region  of  the  Magnetosphere. 

Calculated  Gyroresonant  Energy  at  the  Equator  for 
a  Frequency  of  17.1  kHz,  a  Wave  Normal  Angle  of  60®, 
and  a  Base  Density  of  9200/cm^ . 


Pitch  Angle  Distribution  of  Electrons  in  the  235  keV 
Channel  Observed  in  the  L=1.66  Drift  Shell  Compared 
with  the  Local  Bounce  Loss  Cone  and  the  100  km  Atmospheric 
Loss  Cone  at  45°  East  Longitude . 


Similar  to  Fig.  3,  except  that  the  Electrons  Are 
in  the  L=2.28  Drift  Shell  and  the  Solid  Lines 
Represent  the  Atmospheric  Loss  Cone  at  56° . 


Similar  to  Figs.  3  and  4,  but  for  the  L=4.48  Drift 
Shell . 

Diagram  of  the  Electron  Energies  as  a  Function  of  L, 
Which  Should  Be  Observable  by  the  S3-3  Satellite  for 
the  Pass  of  Fig.  1  under  the  Assumption  that  UMS 
Precipitated  Electrons  of  All  Energies  into  the  Drift 
Loss-Cone  during  Two  Specific  Transmission  Periods..., 


Calculations  of  Resonance  between  Energetic  Electrons 
and  17.1  kHz  Waves . 
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INTRODUCTION 


The  precipitation  of  magnetospheric  electrons  by  gyroreso- 
nant  interactions  with  electromagnetic  waves  from  ground-based 
VLF  transmitters  has  been  a  topic  of  growing  interest  for  almost 
a  decade.  Vampola  [1977a]  presented  observations  of  L-dependent 
precipitating  energetic  electrons  (Eg>100  keV)  in  the  slot  region 
of  the  magnetosphere  which  were  correlated  in  longitude  with  the 
location  of  powerful  ground-based  VLF  transmitters.  In  an  L-de¬ 
pendent  precipitation  event,  the  higher  energy  electrons  precipi¬ 
tate  at  lower  L  in  an  energy  vs  L  relationship  that  can  be  shown 
to  be  consistent  with  a  single  resonant  wave  frequency.  Vampola 
and  Kuck  [1978]  used  observations  of  inner-zone  electrons  in  the 
drift  loss  cone  with  a  trace-back  technique  to  locate  the  actual 
transmitters  which  were  causing  the  precipitation.  Koons  et  al 
[1961]  used  structures  in  the  inner-zone  drift  loss  cone  electron 
distributions  and  synoptic  VLF  data  to  correlate  electron  precip¬ 
itation  events  with  specific  transmission  sequences  from  the  VLF 
transmitters  UMS  and  NWC.  Imhof  et  al  [1983]  correlated  observa¬ 
tions  of  precipitating  electrons  in  the  bounce  loss  cone  of  the 
outer  edge  of  the  inner  zone  with  transmitting  patterns  from  a 
ground-based  VLF  transmitter.  Vampola  [1983]  used  satellite  ob¬ 
servations  of  the  depletion  of  electron  fluxes  transported  across 
the  slot  region  during  a  magnetic  storm  to  determine  the  lifetime 
of  inner-zone  electrons  against  precipitation  by  VLF  transmis¬ 
sions  in  order  to  determine  the  relative  importance  of  VLF  trans- 
m.issions  on  the  morphology  of  magnetospheric  electrons.  He  con¬ 
cluded  that  the  location  and  shape  of  the  outer  edge  of  the  inner 
zone  are  controlled  by  ground-based  VLF  transmissions.  Vampola 
[1987]  used  high-resolution  pitch  angle  measurements  of  precipi¬ 
tating  electrons  in  the  slot  region  to  determine  the  precipita¬ 
tion  pattern  around  a  VLF  transmitter  and  made  resonance  calcula¬ 
tions  which  indicated  the  interaction  region  for  slot  precipita¬ 
tion  was  well  away  from  the  equator. 

While  it  is  now  generally  accepted  that  man-made  waves  are 
instrumental  in  precipitating  electrons  at  the  outer  edge  of  the 
inner  zone,  there  has  been  much  less  acceptance  of  ground-based 
VLF  transm.itters  being  a  significant  factor  in  the  precipitation 
of  electrons  in  the  slot  region.  Precipitation  of  electrons  in 
the  outer  zone  has  not  been  associated  previously  with  these 
ground-based  VLF  transm.itters. 


In  surveys  made  by  one  of  us  (ALV)  of  large  magnetospher ic 
electron  data  bases  obtained  by  a  number  of  electron  spectrome¬ 
ters  flown  on  various  satellites,  a  recurrent  observation  has 
been  that  electron  precipitation  occurs  in  "events"  rather  than 
in  a  slow  diffusion  process,  especially  in  the  slot  region.  Vam- 
pola  and  Gorney  [1983]  made  a  statistical  study  of  energy  deposi¬ 
tion  in  the  middle  atmosphere  by  precipitating  electrons  and  ob¬ 
served  that  for  even  as  high  as  L=4  the  hemispheric  pattern  of 
precipitation  was  that  which  would  occur  if  electron  scattering 
from  stable  drift  orbits  into  the  loss  cone  were  primarily  large 
angle  scattering  rather  than  a  slow  pitch  angle  diffusion.  While 
it  is  known  that  non-resonant  mechanisms  precipitate  electrons  at 
times  (e.g.,  electrostatic  wave  structures  at  the  outer  edge  of 
the  outer  zone  [Koons  et  al,  1972]),  it  is  appropriate  to  ques¬ 
tion  whether  VLF  transmitters  are  also  instrumental  in  precipi¬ 
tating  electrons  in  the  outer  zone.  In  this  paper,  we  will  show 
that  VLF  transmitters  may,  in  fact,  precipitate  large  numbers  of 
electrons  in  the  outer  zone. 


OBSERVATIONS 


The  data  presented  in  this  study  were  obtained  on  28  July 
1976  by  a  magnetic-focusing  electron  spectrometer  on  the  USAF 
Space  Test  Program  satellite  S3-3  while  the  satellite  was  below 
the  stable  trapping  region  of  the  magnetosphere.  That  is,  any 
electrons  mirroring  at  or  below  the  altitude  of  the  satellite  had 
had  their  pitch-angles  lowered  at  some  location  during  the  time 
they  drifted  from  approximately  30°  east  longitude  (EL)  to  the 
location  of  the  satellite  at  approximately  120°  EL.  West  of  30° 
EL,  the  mirror  altitudes  would  have  been  at  or  below  100  km,  the 
nominal  atmospheric  absorption  altitude.  The  spectrometer  con¬ 
sisted  of  two  separate  analyzing  chambers,  one  covering  the  en¬ 
ergy  range  12  keV  to  162  keV  in  five  differential  energy  channels 
and  another  covering  the  range  235  keV  to  1630  keV  in  seven  dif¬ 
ferential  energy  channels.  The  above  energies  refer  to  the  cen¬ 
troid  of  the  energy  response  for  a  "flat"  (energy-independent) 
spectrum.  The  lower  energy  unit  had  a  collimator  with  a  2°x5o 
half-angle  response,  the  other  had  a  5ox5o  half-angle  collimator. 
Geometric-energy  factors  for  the  lower  energy  chamber  were  sev¬ 
eral  orders  of  magnitude  smaller  than  those  of  the  higher  energy 
cham.ber.  The  direction  of  scan  (using  the  spin-stabilized  space¬ 
craft  as  a  scanning  platform)  was  across  the  2°  direction  (in  the 
plane  of  the  instrument) .  The  combination  of  spin  rate,  nominally 
3  RP.M,  and  read-out  of  accumulators,  16  times  per  second,  pro¬ 
vided  pitch-angle  measurements  at  about  1.1°  intervals.  For  this 
study,  data  from,  the  108,  162,  235,  433,  654,  and  880  keV  chan¬ 
nels  only  are  used.  The  data  from  the  12,  33,  and  70  keV  channels 
were  not  used  because  the  70  keV  detector  was  noisy  and  electrons 
of  12  and  33  keV  may  be  present  on  the  drift  shell  due  to  accel¬ 
eration  of  ionospheric  electrons  up  field  lines  by  electrostatic 
potentials . 

The  event  which  is  utilized  in  this  study  was  characterized 
by  precipitation  peaks  in  the  inner  zone,  the  slot,  and  the  outer 
zone.  Events  of  this  type  are  frequently  seen  whenever  the  satel¬ 
lite  is  in  the  drift  loss  cone  of  the  longitudes  associated  with 
the  transmitter  UMS .  In  this  study,  7%  of  the  orbits  examined  had 
such  com.binations  of  slot  and  outer  zone  precipitation  struc¬ 
tures.  Since  many  of  the  orbits  examined  did  not  include  measure¬ 
ments  in  the  drift  loss  cone,  the  actual  incidence  of  such  pre¬ 
cipitations  is  higher.  The  specific  event  selected  for  analysis 
in  this  study  was  selected  on  the  basis  of:  a)  there  being  a  sta¬ 
tistically  significant  intensity  in  all  three  precipitation  peaks 
in  a  single  channel  (required  for  determining  the  longitude  of 


origin)  ;  b)  the  observation  of  flux  in  two  or  more  channels  for 
each  of  the  peaks  (needed  for  a  drift-rate  analysis  comparison 
with  the  transm.ission  schedule  of  the  transmitter)  ;  and  c)  the 
availability  of  the  transmission  schedule  of  the  suspected  trans¬ 
mitter  for  the  time  of  the  precipitation. 

The  precipitation  event  selected  included  measurable  flux  in 
channels  up  to  the  880  keV  channel,  but  the  235  keV  channel  had 
much  higher  counts  due  to  the  facts  that  the  flux  intensity  de¬ 
creased  sharply  with  energy  above  235  keV  and  that  the  lower  en¬ 
ergy  chamber  had  geometric-energy  factors  much  smaller  than  the 
higher  energy  chamber.  Since  accurate  pitch  angle  determinations 
are  required  for  the  type  of  analysis  which  is  used  here,  the  235 
keV  channel  has  been  selected  for  primary  presentation  and  analy¬ 
sis.  The  other  channels  are  used  to  establish  that  precipitation 
features  discussed  here  were  L-dependent  precipitation  events 
and,  because  drift  rates  are  strongly  dependent  on  electron  en¬ 
ergy,  for  detailed  comparison  with  the  transmitting  sequences  of 
the  presumed  transmitter.  In  these  events,  the  235  keV  channel 
typically  had  count-rates  one  to  two  orders  of  magnitude  greater 
than  any  of  the  other  channels.  The  absolute  upper  and  lower  lim¬ 
its  of  response  of  that  channel  are  320  keV  and  158  keV,  respec¬ 
tively,  with  20%  of  the  response  being  above  265  keV  and  20%  be¬ 
ing  below  200  keV,  both  figures  referring  to  a  "  flat  "spectrum. 
The  geometric-energy  factor  for  the  235  keV  channel  for  a  "flat" 
spectrum  is  6.54  cmi^-ster-keV. 

The  counts/0.0625  sec  from  the  108  keV  through  654  keV  chan¬ 
nels,  converted  to  flux  units  using  the  "flat"  spectrum  geomet¬ 
ric-energy  factors,  are  shown  in  Fig.  1.  A  table  at  the  bottom  of 
the  figure  lists  the  geographic  and  geomagnetic  coordinates  of 
the  satellite  during  the  time  the  data  were  being  obtained.  The 
data  have  not  been  corrected  for  background  effects.  The  dashed 
lines  in  the  figure  are  the  fluxes  measured  at  the  same  L  values 
at  higher  altitudes  (between  3600  km  and  7000  km,  the  higher 
value  at  higher  L)  in  the  opposite  hemisphere  and  are  included 
for  comparative  purposes. 

Note  that  the  data  of  Fig.  1  cover  inner  zone,  slot  and 
outer  zone  regions,  all  in  the  drift  loss  cone.  Any  electrons  ob¬ 
served  are  only  quasi-trapped,  have  been  scattered  down  the  field 
line  somewhere  between  30°  EL  and  their  present  longitude,  and 
will  be  lost  into  the  atmosphere  prior  to  or  near  the  South  At¬ 
lantic  Anom.aly.  The  numbers  3  through  5  in  the  235  keV  panel  re¬ 
fer  to  the  locations  at  which  data  were  obtained  for  presentation 
in  Figs.  3  through  5. 


Figure  1.  Electron  fluxes  from  the  S3-3  108  through  654  KeV  channels  obtained  on  28  July  1976  while 
the  satellite  was  below  the  stable  trapping  region  of  the  magnetosphere.  All  electrons  ob¬ 
served  during  this  pass  will  be  lost  into  the  atmosphere  before  completing  one  drift  orbit. 
Any  electrons  observed  during  this  pass  must  have  entered  the  drift  loss  cone  east  of 
about  30=  east  longitude.  The  dashed  lines  show  stably  trapped  fluxes  measured  in  the 
opposite  hemisphere  at  altitudes  between  3600  km  (L=1.5)  and  7000  km  (L=7.0).  TTie 
numbers  3-5  refer  to  the  regions  of  data  plotted  in  Figs.  3-5. 
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Figure  2  is  a  plot  of  the  location  in  L  of  the  peak  precipi¬ 
tation  observed  in  the  inner  zone  in  the  108,  162,  235,  and  433 
keV  channels  and  in  the  slot  at  108,  162,  235,  433,  and  654  keV 
plotted  at  their  nominal  energies.  The  bars  indicate  the  total 
energy  span  of  each  channel.  The  solid  line  is  a  calculation  of 
the  nominal  equatorial  resonant  energy  vs  L  for  a  wave  of  fre¬ 
quency  17.1  kHz  (the  UMS  frequency)  with  a  wave-normal  angle  of 
6C^,  a  base  density  of  9200/cm3,  and  using  a  simple  exponential 
decrease  of  electron  density  with  altitude  in  a  dipole  field.  The 
wave-normal  angle  was  selected  on  the  basis  of  a  previous  ray¬ 
tracing  analysis  in  the  inner  zone  for  UMS  [Koons  et  al,  1981]. 
The  base  density  was  selected  as  being  the  best  fit  to  the  ray¬ 
tracing  analysis  for  the  235  keV  channel  data  which  is  presented 
later  in  the  Resonance  Calculations  section.  There  is  no  question 
that  both  the  inner-zone  and  slot  precipitation  are  L-dependent . 
The  fit  to  the  inner-zone  precipitation  peaks  is  not  exact  be¬ 
cause  a  relatively  crude  resonance  calculation  was  performed  for 
Fig.  2  and,  for  reasons  which  will  become  clear  later  in  the  dis¬ 
cussion,  the  energy  of  the  particles  detected  by  a  given  channel 
m.ay  net  be  well  represented  by  using  the  centroid  of  response  of 
that  channel . 

The  first  peak,  labeled  "3"  in  Fig.  1,  is  an  L-dependent 
precipitation  of  the  type  analyzed  in  detail  by  Koons  et  al 
[1961],  Sufficie.nt  flux  for  L-determination  was  observed  in  the 
iC5,  162,  235  and  433  keV  channels,  but  due  to  the  small  geomet¬ 
ric-energy  factors  of  the  lower  channels  only  the  235  keV  channel 
had  a  sufficient  count-rate  for  pitch-angle  determination  with 
accuracy  of  the  type  required  for  use  in  a  traceback-to-longi- 
tude-cf-crigin  technique  [Luhm.ann  and  Vampola,  1977].  Analysis  of 
a  large  numlcer  cf  similar  events  by  Vam.pola  and  Kuck  [1978]  indi¬ 
cated  that  inner- zone  precipitation  with  energies  in  excess  of 
4GC'  keV  IS  rare,  so  it  is  probably  not  significant  that  few  elec¬ 
trons  were  observed  at  energies  above  the  235  keV  cha.nnel. 

The  traceback-tc-longitude-of-origin  technique  is  applicable 
to  particles  observed  within  the  drift  loss  cone  under  the  as- 
sumipticn  that  the  particle  pitch  angle  distribution  was  produced 
by  1  nteract i o.n  with  the  atmosphere  at  the  longitude  of  precipita¬ 
tion  and  has  not  been  substantially  miodified  during  the  time  the 
particles  have  drifted  from  their  longitude  of  precipitation  to 
the  longitude  of  observation.  The  procedure  requires  transforming 
the  local  pitch  angle  to  equatorial  pitch  angle  and  then  matching 
the  equatorial  pitch  angle  cutoff  to  the  equatorial  atmospheric 
loss  cone  angle  at  som.e  previous  longitude,  with  the  requirem.ent 
that  no  intervening  longitude  has  an  atmospheric  loss  cone  angle 
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Figure  2,  Calculated  gyroresonant  energy  at  the  equator  for  a  frequency  of  17.1  kHz,  a  wave  normal  an¬ 
gle  of  60-,  and  a  base  density  of  92(X)/cm3.  The  dots  are  the  inner  zone  precipitation 
peaks  from  the  108,  162,  235,  and  433  keV  channels  and  the  slot  precipitation  p>eaks  at 
108,  162,  235,  433,  and  654  keV  plotted  at  their  nominal  energy.  The  bars  indicate  the 
total  energy  span  of  the  channels. 


Figure  3  presents  the  raw  counts  for  one-half  spin  period  of 
the  satellite,  normalized  to  the  85°-95o  pitch  angle  average, 
plotted  as  a  function  of  pitch  angle  for  the  data  of  peak  3  in 
Fig.  1.  The  L-value  of  this  data,  1.66,  is  similar  to  those  of 
the  events  reported  by  Vampola  and  Kuck  [1978],  Koons  et  al 
[1981],  and  by  Imhof  et  al  [1981]  in  which  electrons  in  this  en¬ 
ergy  range  were  being  precipitated  in  the  inner  zone  by  VLF 
transmitters.  The  smaller  intensity  in  this  peak  compared  to 
peaks  4  and  5  is  probably  a  matter  of  instrumental  response,  for 
reasons  to  be  discussed  later,  and  should  not  be  taken  as  an  in¬ 
dication  that  fewer  electrons  were  precipitated  in  the  inner 
zone.  The  dashed  lines  show  the  response  of  the  instrument  to  a 
distribution  in  which  the  flux  is  zero  within  the  local  atmo¬ 
spheric  loss  cone  and  isotropic  beyond  it.  The  observed  particle 
pitch  angle  distribution  is  much  narrower  than  the  local  loss 
cone,  indicating  that  the  pitch  angle  distribution  is  not  deter¬ 
mined  by  local  conditions  but  rather  by  scattering  and  atmo¬ 
spheric  cutoff  conditions  at  some  previous  longitude  (s)  •  The 
solid  lines  represent  the  instrument  response  to  a  distribution 
defined  by  isotropic  precipitation  at  45°  EL,  the  western  end  of 
the  UMS  precipitation  region  [Koons  et  al,  1981].  The  data  are 
fit  quite  well  by  the  assumption  of  a  45°  EL  origin.  For  a  better 
fit  to  the  slope  of  the  pitch  angle  distribution,  one  would  have 
to  take  into  account  the  entire  longitudinal  precipitation  pat¬ 
tern  that  is  contributing  to  the  data  points,  as  was  done  in  Vam¬ 
pola  [1987].  The  magnetic  field  model  used  for  tracing  purposes 
for  this  and  the  other  peaks  was  the  DGRF  1975  [lAGA,  1981], 
Epoch  1976,  internal  field  model.  No  external  field  model  was 
used.  It  can  be  concluded  from  the  L-dependent  nature  of  the  pre¬ 
cipitation,  the  fit  of  the  L  vs  Energy  dependency  to  the  reso¬ 
nance  calculation,  and  from  the  longitude  of  its  origin  that  the 
flux  in  peak  3  was  undoubtedly  precipitated  by  the  ground-based 
VLF  transmitter  UMS.  The  L-value  of  UMS  varies  from  1.73  to  1.77, 
depending  on  whether  44°  EL  or  56°  EL  and  whether  0  or  100  km,  al¬ 
titude  is  used  for  the  determination.  Calculations  presented  in 
the  Resonance  Calculations  section  were  performed  with  the  UMS 
frequency,  17.1  kHz,  and  location,  44°  EL,  46.2°  N  [Watt,  1967]. 

Since  the  particles  trace  back  to  the  western  end  of  the  UMS 
precipitation  region,  it  appears  that  the  particles  may  be  the 
tail  end  of  a  longitudinal  distribution  due  to  a  discrete  trans- 
miission  period.  When  the  transmitter  turns  off,  precipitation 
stops  and  the  entire  longitudinal  pattern  of  precipitation  drifts 
eastward,  with  different  drift  rates  for  different  energies  and 
L-shells.  (The  drift  rate  is  relatively  insensitive  to  pitch  an¬ 
gle  for  the  small  variations  in  pitch  angle  of  concern  here.)  We 
shall  address  this  topic  again  in  the  Analysis  section. 
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Figure  3.  Pitch  angle  distribution  of  electrons  in  the  235  KcV  channel  observed  in  the  L=1.66  drift  shell 
compared  with  the  local  bounce  loss  cone  (dashed  lines)  and  the  100  km  atmospheric  loss 
cone  at  45=  east  longitude  (solid  lines).  The  pitch-angle  distributions  match  the  100  km 
atmospheric  loss  cone  at  the  westernmost  longitude  of  UMS  precipitation.  See  text  for 
details. 
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Figure  4  is  sir.ilar  to  Fig.  3  but  is  slot  data  at  a  higher 
L-value  than  the  in.ner-zcne  precipitation.  The  L-value  is  similar 
to  that  of  the  observations  of  Imhof  et  al  [1983],  but  the  energy 
is  significantly  higher.  The  observations  of  Imhof  et  al  [1983] 
wore  at  quite  low  energy  and  were  the  low  energy  extension  (to 
higher  L)  of  the  inner-zone  equatorial  region  electron  precipita¬ 
tion  by  ground-based  VLF  transmitters  (Vam.pola  and  Kuck  [1978]). 
The  peak  of  Fig.  4  was  again  an  L-dependent  peak  as  can  be  seen 
from.  Fig.  2,  with  strong  fluxes  (400  counts/sec)  in  the  235  keV 
channel,  significant  fluxes  in  the  433  keV  channel,  and  traces  of 
flux  in  the  108,  162,  654  and  880  keV  channels.  The  dashed  lines 
are  again  as  defined  in  Fig.  3,  but  in  Fig.  4  the  solid  lines 
represent  the  atmospheric  loss  cone  at  56°,  the  eastern  end  of 
the  UMS  precipitation  region.  These  match  the  data  points,  giving 
an  unam.biguous  fit  to  the  eastern  end  of  the  UMS  precipitation 
region.  Since  the  eastern  end  of  the  precipitation  region  has  the 
sm.allest  equatorial  pitch  angle  cutoff,  the  simple  comparison 
with  pitch  angle  cutoffs  being  made  here  should  always  trace  back 
to  the  easter.nm.ost  longitude  of  the  precipitation  region  that  is 
currently  represented  in  the  drifting  particles  at  the  longitude 
of  the  satellite.  Figures  3  and  4  indicate  that  a  single  trans- 
m.itter,  U”S ,  was  probably  precipitating  electrons  simultaneously 
on  both  an  inner-zone  field  line  and  a  slot  field  line.  We  will 
s.how  later,  using  drift-rate  analyses  and  the  operating  times  of 
UMS,  that  this  was  indeed  the  case. 

Figure  5  presents  the  L=4.48  data  (peak  5  of  Fig.  1)  in  a 
fcrm.at  sim.ilar  to  Fig.  3.  For  this  event,  flux  was  observable  in 
the  235,  433  and  654  keV  channels.  These  particles  were  precipi¬ 
tated  east  of  30°  EL  within  the  previous  18  minutes  (the  time  for 
electrons  with  this  energy,  pitch-angle,  and  L- she 11  to  drift 
from.  30°  EL,  the  westernmost  location  on  this  drift  shell  which 
can  be  observed  by  the  satellite  at  this  altitude  and  longitude) . 
Due  to  the  broad  region  over  which  the  235  keV  channel  is  detect¬ 
ing  electrons,  it  is  difficult  to  assign  a  single  L-value  for  the 
peak.  The  match  between  the  solid  lines  and  the  data  points  in 
Fig.  5  again  indicates  that  the  electrons  trace  back  to  the  east¬ 
ern  end  of  the  UMS  precipitation  region  and  were  precipitated  in 
the  vicinity  of  U.MS .  For  both  the  L=2.28  and  L=^.48  precipitation 
events,  the  trace-back  is  accurate  to  within  several  degrees  of 
longitude  of  the  eastern  end  of  the  UMS  precipitation  region.  Two 
degrees  in  pitch  angle  are  equivalent  to  about  one  degree  in  lon¬ 
gitude  around  56°  EL  for  this  geometry.  The  fact  that  this  pitch 
angle  distribution  maps  back  to  the  location  of  UMS  at  a  time 
when  UMS  was  precipitating  electrons  from,  the  inner  zone  and  from, 
the  slot  is  indicative  of  these  ou^er  zone  electrons  being  pre- 
C-cpitated  by  UMS  even  though  there  is  no  obvious  L-dependency . 


Figure  4.  Similar  to  Fig.  3,  except  that  the  electrons  are  in  the  L=2.28  drift  shell  and  the  solid  lines  rep¬ 
resent  the  atmospheric  loss  cone  at  56°.  The  fit  to  the  atmospheric  loss  cone  at  the  east¬ 
ernmost  longitude  of  UMS  precipitation  (56°)  is  unambiguous. 
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Figure  5.  Similar  to  Figs.  3  and  4,  but  for  the  L=4.48  drift  shell.  Again,  the  data  fit  the  atmospheric  loss 
cone  at  the  longitude  of  UMS. 
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ANALYSIS 


Because  of  the  highly  surprising  nature  of  the  conclusion, 
above,  that  UMS  was  precipitating  electrons  from  the  outer  zone, 
additional  supportive  evidence  is  necessary,  even  though  the  di¬ 
rect  evidence  is  strong.  Direct  comparison  of  the  transmitting 
schedule  of  UMS  with  the  particle  observations  would  be  m.ost 
helpful.  If  UMS  were  not  operating,  or  if  its  operating  schedule 
were  not  consistent  with  the  particle  observations,  another  ex¬ 
planation  would  have  to  be  sought  for  the  precipitations.  Fortu¬ 
nately,  the  transmission  sequence  of  UMS  at  the  appropriate  time 
is  available  for  comparison  with  the  drift  rates  of  the  elec¬ 
trons.  One  can  calculate  the  drift  rate  for  electrons  of  various 
energies  on  the  various  L-shells  using  the  relationship  [Schulz 
and  Lanzerotti,  1974]: 


2K/CI2  =  -{3L/27rY)  (y2-i)  (c/a)  2  (mg/gBo)  D  (y) /T  (y) 

where  ^2  the  drift  frequency,  y  is  the  relativistic  mass  ratio 
of  the  electron,  c  is  the  velocity  of  light,  a  is  the  earth  ra¬ 
dius,  L  is  in  units  of  earth-radius,  mg  is  the  mass  of  the  elec¬ 
tron,  Bq  is  the  value  of  the  earth's  magnetic  field  at  the  sur¬ 
face  of  the  earth  at  the  equator,  q  is  the  charge  on  the  elec¬ 
tron,  y  is  sina  where  a  is  the  equatorial  pitch  angle,  and  D(y) 
and  T(y)  are  given  by: 


D(y)  =  Vi2{5.5208-.4381y-.6397  (y  Iny  +  y-'^)} 

T  (y)  =1 . 38  02-.31985  (y  +  y-'^) 

By  making  use  of  the  fact  that  electrons  with  various  ener¬ 
gies  on  a  given  L-shell  drift  with  different  rates,  the  data  from 
each  energy  channel  become  a  separate  test  of  the  validity  of  the 
assignment  of  UMS  as  the  cause  of  the  precipitation.  Each  opera¬ 
tional  period  of  UMS,  precipitating  all  energies  at  the  same  time 
at  the  same  longitudes,  produces  a  series  of  longitudinally- 
drifting  electron  bunches  characterized  by  a  finite  longitudinal 
length  (determined  by  the  length  of  time  the  transmitter  is  radi¬ 
ating  and  the  longitudinal  width  of  the  interacting  region) 
drifting  with  a  rate  that  is  energy  dependent.  At  a  given  later 
time,  different  energies  will  be  found  at  different  longitudes. 
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For  the  comparison  with  the  UMS  transmission  schedule,  data 
were  extracted  by  Rayspan  analysis  from  synoptic  VLF  recordings 
miade  at  the  time  the  satellite  data  were  obtained.  The  receiver 
consisted  of  a  vertical  loop  antenna  of  100  square  meters  area 
with  the  plane  of  the  loop  oriented  at  70°  east  of  north,  feeding 
a  preamplifier  covering  the  frequency  range  from  0.5  to  30  kHz. 
The  preamplifier  output  was  telemetered  10  km  to  the  University 
of  Otago  campus  by  UHF  link  for  analog  recording  on  magnetic 
tape.  Timing  information  accurate  to  a  few  milliseconds  was  pro¬ 
vided  by  NASA-36  timecode  recorded  on  a  separate  track. 

In  order  to  ascertain  when  UMS  was  transmitting,  dynamic 
spectra  from  this  tape  in  the  range  from  10-20  kHz  were  recorded 
on  film  using  a  modernized  Rayspan  real-time  spectrum  analyzer. 
This  gave  a  detection  limit  for  signals  arriving  from  the  direc¬ 
tion  of  UMS  of  about  60  microvolts  per  meter  while  the  UMS  trans- 
m.ission  was  determined  to  have  a  signal  strength  of  about  200  mi¬ 
crovolts  per  meter.  This  indicates  a  radiated  power  level  of  sev¬ 
eral  hundred  kW  based  on  sub-ionospheric  propagation  calculations 
following  Watt  [1967],  This  power  level  is  in  agreement  with  pub¬ 
lished  figures  (315  kW,  Watt  [1967])  for  UMS. 

Using  UMS  transm.ission  times  from  the  Rayspan  analysis  and 
assur.ing  that  all  energies  of  electrons  are  precipitated  by  UMS 
at  all  L-values,  Fig.  6  diagrams  the  electron  energies,  as  a 
function  of  L,  that  should  be  observable  by  the  S3-3  as  it 
crosses  the  various  L-shells  in  the  pass  plotted  in  Fig.  1.  The 
first  transmission,  which  ended  at  12:59:  O’!  UT,  is  diagramed  with 
both  a  56°  and  a  45°  cutoff,  since  the  data  of  Fig.  3  indicate 
the  inner-zone  particles  were  precipitated  near  45°.  Electrons 
precipitated ' by  the  first  transmission  with  energies  above  those 
given  by  the  45°  line  would  have  already  drifted  past  the  longi¬ 
tude  of  the  satellite  by  the  time  the  satellite  arrived.  Simi¬ 
larly,  electrons  precipitated  at  56°  with  energies  above  the  line 
m.arked  56°  would  also  have  drifted  past  the  longitude  of  the 
satellite.  For  the  second  transmission,  which  began  at  13:16:09 
UT,  electrons  with  energies  below  the  line  marked  56°  would  not 
yet  have  arrived  at  the  longitude  of  the  satellite.  The  vertical 
dashed  lines  marked  3,  4,  and  5  indicate  the  L-values  of  the  data 
of  F igs .  3-5  . 

We  will  now  compare  the  predictions  of  observable  electron 
energies  as  a  function  of  L  in  Fig.  6  with  the  actual  observa¬ 
tions.  Figure  6  predicts  that  electrons  with  energies  between 
about  12C  keV  and  155  keV  precipitated  at  56°  EL  by  the  first 
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Figure  6.  Diagram  of  the  electron  energies  as  a  function  of  L  which  should  be  observable  by  the  S3-3 
satellite  for  the  pass  of  Fig.  1  under  the  assumption  that  UMS  precipitated  electrons  of  all 
energies  into  the  drift  loss-cone  during  two  specific  transmission  periods.  Lines  marked 
45 :  and  56-  indicate  the  longitude  at  which  the  precipitation  is  assumed  to  have  occurred. 
Dashed  lines  marked  3-5  indicate  the  L- values  of  the  peaks  marked  3-5  in  Fig.  1. 
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trar.srr.issicr.  ar.d  electrons  with  energies  above  about  335  keV  pre¬ 
cipitated  by  the  second  transmission  should  be  observable  at  an  L 
of  1.66.  Additionally,  electrons  with  energies  between  155  keV 
and  195  keV  precipitated  at  increasingly  western  longitudes  to 
45"  should  also  be  observable.  Figure  1  shows  that  electrons  are 
actually  observed  in  the  108,  162,  235,  and  433  keV  channels.  The 
1C8  keV  channel  responds  to  particles  with  energies  up  to  118  kev 
and  trails  off  to  zero  response  at  121  keV.  The  162  keV  channel 
can  respond  to  electrons  with  energies  as  low  as  151  keV.  Thus, 
electrons  should  be  seen  in  the  108  and  162  keV  channels,  but 
with  a  very  reduced  intensity  since  the  electron  energies  are 
near  the  limits  of  the  response  of  the  channels.  This  is  what  is 
seen.  The  235  keV  channel  can  respond  to  the  195  keV  electrons 
precipitated  by  the  first  transmission,  but  again  with  very  re¬ 
duced  intensity  since  the  major  response  of  that  channel  (80%)  is 
above  20C  keV.  For  the  inner-zone  event,  the  2  35  keV  channel 
traces  back  to  45”  EL  instead  of  56=  EL  because  electrons  precip¬ 
itated  farther  east  which  are  in  the  m.ajor  energy  response  of  the 
channel  have  already  drifted  past  the  location  of  the  satellite. 
The  chan.nel  is  even  less  sensitive  to  electro.ns  with  energies 
sufficiently  low  to  have  been  precipitated  east  of  45"  and  still 
be  intercepted  by  the  satellite. 

The  second  transmission  might  also  have  contributed  to  the 
response  of  the  instrument  at  L=1.66.  The  235  keV  channel  has  a 
very  s.mall  response  at  energies  up  to  4  30  keV,  with  perhaps  1%  of 
its  total  response  being  above  365  keV.  The  very  slight  response 
in  the  433  keV  channel  (which  could  respond  only  to  electrons 
recipitated  by  the  second  transmission)  and  the  lack  of  a  re- 
pcnse  in  the  higher  energy  channels  is  undoubtedly  due  to  the 
act  that  significant  numbers  of  electrons  with  energy  above  400 
eV  are  not  precipitated  in  the  inner  zone  by  UMS  [Vamipola  and 
uck,  19^6:.  Cn  the  basis  of  the  transmitting  schedule  of  UMS, 
lectrcn  drift  tim.es,  and  inner-zone  precipitation  energy  spec- 
ra,  electrons  should  be  seen  in  t.he  108,  162,  235,  and  possibly 

he  433  keV  channels.  That  is  precisely  what  is  observed.  The 
rift-tim.e  analysis  also  properly  explains  t.he  traceback  to  45^ 
L  for  peak  3. 

The  slot  precipitation  analysis  is  similar.  For  this  event, 
electrons  were  observable  in  all  channels  up  to  the  654  keV  chan¬ 
nel.  Electrons  precipitated  by  the  first  transmission  should  be 
observable  in  the  105  and  162  keV  channels,  while  electrons  above 
3C0  keV  from,  the  second  transmission  should  be  seen  in  the  other 
channels.  The  2  35  keV  cha.nnel  has  20%  of  its  response  above  320 
ke'.'  and  was  undoubtedly  m.easuring  electro.ns  above  30C  keV,  al¬ 
though  t.ne  m;easurea  intensity  seem.s  high  for  a  channel  whose  m.a- 
;or  response  ^s  below  this  energy.  The  higher  energy  channels 
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were  seeing  particles  of  appropriately  higher  energy  which  were 
precipitated  by  the  second  transmission.  The  drift  analysis  based 
on  the  known  transmission  times  and  location  again  properly  ac¬ 
counts  for  details  of  the  particle  observations.  This  analysis, 
together  with  the  fact  that  UMS  was  currently  precipitating  in¬ 
ner-zone  electrons,  adds  strong  evidence  to  previous  claims  that 
ground-based  VLF  transmitters  are  responsible  for  significant 
precipitation  in  the  slot  [Vampola,  1977a;  Vampola,  1987].  An  im¬ 
portant  point  to  consider  is  that,  unlike  inner-zone  VLF  trans¬ 
mitter-induced  precipitation,  slot  precipitation  by  ground-based 
VLF  transmitters  precipitates  electrons  with  energies  well  above 
400  keV. 

Ground-based  VLF  transmitter  precipitation  of  electrons  from 
the  inner  zone  is  now  generally  acknowledged  by  the  scientific 
comumunity;  slot  precipitation  produced  by  the  same  transmitters 
has  received  much  less  acceptance.  Precipitation  of  outer  zone 
electrons  by  these  same  transmitters  has  not  been  considered  pre¬ 
viously.  If  such  a  mechanism  is  operating,  the  implications  of  it 
to  the  morphology  of  the  electron  population  of  the  magneto¬ 
sphere,  and  our  understanding  thereof,  are  profound.  Following 
the  analyses  used  for  the  inner  zone  and  slot  observations,  we 
can  similarly  analyze  the  outer  zone  precipitatio.ns  seen  in  Fig. 
1 . 

Figure  6  shows  that  even  the  lowest  energy  electrons  from 
the  first  transmission  would  have  drifted  past  the  longitude  of 
the  satellite  prior  to  the  arrival  of  the  satellite  at  the  longi¬ 
tude  of  observation.  Only  electrons  with  energies  above  210  keV 
precipitated  by  the  second  transmission  would  be  observable  by 
the  satellite  at  the  time  and  longitude  of  observation.  Electrons 
were  actually  seen  in  the  235,  433,  and  654  keV  channels.  A  very 
important  point  is  that  no  electrons  were  observed  in  the  108  and 
162  keV  cha.tnels.  Electrostatic  wave  scattering  [Koons  et  al, 
1972]  or  magnetospher ic  hiss  and  chorus  would  not  be  energy  se¬ 
lective  and  would,  in  fact,  be  expected  to  scatter  low  energy 
electrons  m.cre  efficiently  than  high  energy  electrons.  The  high 
altitude  measuremients  in  the  outer  zone  on  the  same  orbit  showed 
electrons  with  energies  in  the  108  and  162  keV  channel  range  were 
present  in  significant  numbers.  In  the  drift  loss  cone,  electrons 
are  observed  specifically  in  the  235,  433,  and  654  keV  channels 
not  because  only  those  energies  are  being  precipitated,  but  be¬ 
cause  lower  energy  electrons  from  the  first  transmission  had  al¬ 
ready  drifted  past  the  longitude  of  observation  and  those  from 
the  seccnd  transmission  have  not  yet  had  time  to  drift  to  the 
longitude  of  observation. 
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We  have,  in  this  event,  very  significant  precipitation  of 
energetic  electrons  in  the  outer  zone  which  is  completely  ac¬ 
counted  for  in  every  detail  by  using,  as  its  source,  precipita¬ 
tion  over  UMS  while  UMS  was  transmitting.  Conversely,  the  energy 
spectrumi  of  the  precipitation  is  unlike  that  which  would  be  pro¬ 
duced  by  any  other  known  source.  The  observation  was  obtained  in 
close  association  with  observations  of  precipitation  by  that  same 
transm.itter  in  the  inner  zone  and  in  the  slot.  This  and  the 
agreement  between  the  transmission  schedule  and  the  drift-rate 
analysis,  we  believe,  unambiguously  identifies  this  outer  zone 
precipitation  as  being  caused  by  UMS  and  establishes  a  previously 
unknown  phenomenon:  precipitation  of  energetic  outer  zone  elec¬ 
trons  by  ground-based  VLF  transmitters.  Additional  support  for 
the  conclusion  that  UMS  produces  energetic  electron  precipitation 
in  the  outer  zone  can  be  gained  from  resonance  calculations. 
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In  addition  tc  establishing  that  resonance  can  occur  between 
'.i  St  ler-rr.ode  VLF  waves  with  the  frequency  cf  I'MS  and  gyrating 
.ectrons  of  the  energy  observed  in  this  study,  another  rationale 
)r  performing  resonance  calculations  is  the  determination  of  the 
robable  location  on  the  field  line  for  the  gyroresonant  interac- 
-on.  A  numb)er  of  investigations  of  wave-particle  interactions  in 
le  magnetosphere  have  been  made  previously  with  various  assum.p- 
■ons.  For  instance,  Helliwell  et  al  [1975]  assumed  interactions 
-th  ducted  waves  in  the  equatorial  region  because  such  condi- 
-cns  allowed  a  m.ore  tractable  calculation.  Inan  et  al  [1978] 
.mulated  interactions  between  coherent  waves  and  low  energy 
.ectrons  and  found  that  the  interaction  is  linear  when  the  wave 
'.tensity  is  low  com.pared  to  the  transverse  energy  of  the  parti- 
.e  (3  pT  for  1  keV  electrons  at  10°  pitch-angle)  or  for  interac- 
-or.s  away  frcm>  the  equator  where  the  inhomiogene ity  of  the  medium, 
jm.inates .  boons  et  al  [1961]  utilized  wave-tracing  of  ncn-ducted 
ive  propagation  to  an  inner-zone  interaction  region  near  the 
ruatcr  in  com.parisons  with  L-dependent  energetic  electron  pre- 
-pitations.  Vampcla  [1977a]  hypothesized  wave-particle  interac- 
■  ons  in  the  slot  region  m.uch  lower  on  the  field  line  but  did  not 
resent  detailed  calculations. 


f  calculations  have  been  consistent  with  cbserva- 
:  al  [1961]  calculated  that  a  field  intensity  of  3 
:nt  to  produce  the  scattering  they  observed,  while 
cf  electron  interaction  with  coherent  waves  re- 
et  ai  [1963]  indicated  that  detectable  electron 
-Id  be  accom^pl i shed  with  as  little  as  1  pT.  They 
energy,  a  few  keV,  electrons  and  assum.ed  ducting. 

reported  in-situ  m.ea su re.men t s  near  the  equator 
cT  levels  from,  the  Siple  transm.itter  before  magne- 
ificatic.n  and  also  reported  observat  ion  of  an  am- 
ricr  to  its  first  passage  through  the  equatorial 
ting  that  wave-particle  interactions  could  occur 
the  eouator. 


Eecause  the  ionosphere  is  a  miagnetized  plasm.a,  its  proper¬ 
ties  are  anisotropic  for  radio  waves  reflecting  off  of  it  or 
rotagating  through  it.  Furthermiore,  it  is  generally  highly  at- 
enuating.  As  a  result,  realistic  calculations  of  wave  propaga- 
lon  have  not  been  m.ade  in  analytic  form,.  Instead,  ray-tracing 
techniques  are  used.  The  advantage  of  ray-tracino  is  that  one  can 
describe  the  propagation  medium,  in  either  continuous  or  discrete 
form  and  make  relatively  good  predictions  of  the  resultant  wave 
r  r  op  aoat  1  o.".  vector  and  am.plitude  without  resortina  to  ducts.  Con- 


() 


stuaies  nave  sr.cwn  wave  propagation  between  ir.id-latitude 
te  points.  buo:ting  is  usually  given  as  the  explanation, 
propagation  requires  that  an  otherwise  reflecting  iono- 
te  sufficiently  tra.nsparent  to  waves  for  significant  en- 
leak  throuah  to  the  duct. 
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The  georr.etry  of  the  calculation  of  resonance  between 
ler-rr.ode  waves  and  energetic  electrons  is  as  follows:  VbF 
from,  a  ground-based  transmitter  in  the  northern  hemisphere, 
gating  upward  and  southward  in  an  unducted  mode,  interact 
energetic  electrons  travelling  northward  and  dow.nward  along 
ocal  magnetic  guiding  center.  Some  of  the  electrons  interact 
esonantly  with  the  waves,  giving  up  some  of  their  transverse 
y.  I.n  the  potential  well  of  the  wave,  some  transverse  energy 
so  converted  to  parallel  energy  [Inan  et  al,  1983].  The  wave 


and  the  local  pitch  angle  of  the  particle  is  re- 


If  sufficient  energy  transfer  occurs  or  if  sufficient  en- 
s  lost,  the  pitch  angle  may  be  reduced  sufficiently  to 
the  particle's  mirror  altitude  to  below  the  atmospheric 
(ICC  km) .  Particles  with  lowered  pitch  angles  which  escape 
lest  into  the  atmosphere  at  the  longitude  of  the  transmit- 
ift  eastward  where  they  may  be  observed  by  the  S3-3  instru- 
no  are  ultim.ately  lost  into  the  atmosphere  at  a  longitude 


:he  particle. 


resonance  condition  for  an  unducted  whistler-mode  wave 
ncy  OJ  and  wave  normal  0  interacting  with  a  relativistic 


a  local  pitch  angle  a  is  given  by 


CO  ( 1- ppcosBcosa)  =  mQ  (1- P") 


ne  electron  velocity,  is  the  cold  elec- 
is  the  order  of  the  resonance,  and  [i  is  the 
iven  bv 


p  =  {  1  -  COj^-ZlCO  (il. 


)S0  -  CO)]}-' 
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where  cOp^  is  the  plasrr.a  frequency,  squared,  and  is  given  by 


C0p2  =  4Kqg2Ne  [  (l/mg)  +  (1/mp)  ] 

where  Ng  is  the  local  plasma  density,  qg  is  the  electron  charge, 
and  is  the  rest  mass  of  the  electrons/ions  making  up  the 
plasma.  In  a  realistic  plasmaspheric  model,  a  combination  of  pro¬ 
tons  and  singly-ionized  oxygen  ions  is  usually  used  and  an  appro¬ 
priately  weight-.;d  mass  is  used  for  computing  the  plasma  fre¬ 
quency  . 

Figure  7  shows  the  resonant  energy  as  a  function  of  altitude 
for  first  order  resonance  with  unducted  waves  for  the  three  re¬ 
gions  of  precipitation  discussed  in  this  study.  The  frequency  was 
selected  to  correspond  to  that  of  UMS;  the  base  density  for  a 
diffusive  equilibrium  model  of  the  plasmasphere  was  selected  to 
give  a  best  fit  to  the  inner-zone  precipitation  data;  and,  the 
ratio  was  set  at  a  nominal  90%/10%  at  a  temperature  of  1600 
°K  (following  Koons  et  al  [1981])  .  The  L  values  selected  are 
those  which  are  representative  of  the  center  of  the  precipitation 
peaks  in  the  235  keV  channel.  The  index  of  refraction  and  wave 
normal  angle  were  determined  using  a  wave-tracing  procedure  based 
on  the  work  of  Burtis  [1973]  as  modified  by  Edgar  [1976]  and 
Thomson  [1976].  The  Thomson  modification  is  a  provision  for  vary¬ 
ing  the  latitude  at  which  the  light  ion  and  heavy  ion  densities 
are  equal.  The  modification,  called  the  Transition  Level  Gradient 
(TLG)  m.odel,  provided  a  better  fit  to  the  data  of  ISIS  II  than 
the  original  Burtis  codes  as  modified  by  Edgar.  The  base  density 
for  a  diffusive-equilibrium  model  of  the  ionosphere/plasmasphere 
was  selected  on  the  basis  of  a  best  fit  to  the  235  keV  channel 
data  at  L=1.66.  The  selected  equatorial  pitch  angle  of  the  elec¬ 
tron  is  that  which  just  permits  a  particle  to  survive  drift 
through  the  South  Atlantic  Anomaly.  This  value  was  selected  on 
the  basis  that  it  requires  the  minimal  change  in  pitch  angle  (and 
energy  loss)  to  precipitate  it  into  the  bounce  or  drift  loss 
cone . 

The  resonance  calculation  is  not  completely  consistent  in¬ 
ternally,  since  the  ray-tracing  codes  used  a  dipole  field  model 
while  the  resonance  calculation  followed  the  field  line  using  the 
sarnie  field  model  as  was  used  in  the  drift-rate  and  longitude-of- 
origin  studies:  DGRF  1975  [lAGA,  1981],  Epoch  1976.  The  ray  trac¬ 
ing  procedure  assumed  that  waves  were  injected  vertically,  giving 
an  initial  wave-norm.al  angle  equal  to  the  local  dip-angle  of  the 
field  line.  The  waves  did  not  follow  the  field  line,  but  a  series 
of  traces  at  slightly  different  latitudes  (plus  or  minus  a  few 
degrees  around  the  latitude  of  the  transmiitter )  were  used  to  de¬ 
termine  the  wave-norm.al  angle  and  index  of  refraction  at  a  numi^er 
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Figure  7.  Calculations  of  resonance  between  energetic  electrons  and  17.1  kHz  waves.  Resonance  regions 
are  show'n  which  compare  with  the  data  of  Figs.  3,  4,  and  5.  The  L=1.66  resonance  for 
195  KeV  electrons  is  at  the  equator.  The  other  resonance  locations  are  well  off  the 
equator  at  altitudes  of  5000  to  6000  km. 
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cf  points  along  the  field  line.  In  the  resonance  calculation,  in¬ 
terpolation  was  then  used  to  determine  the  local  wave-normal  and 
local  index  of  refraction. 

Figure  7  shows  that  first  order  resonance  conditions  for 
electrons  of  the  energy  (-300  KeV)  and  L-value  (2.28)  of  Fig.  4 
are  obtained  in  the  vicinity  of  5000  km,  well  away  from  the  equa¬ 
tor.  For  resonance  to  occur  near  the  equator  with  the  particles 
observed  by  the  S3-3  at  this  L-value,  transmitter  frequencies  be¬ 
low  2C00  Hz  would  be  required,  well  below  operational  VLF  trans¬ 
mitter  frequencies.  For  resonance  to  occur  at  the  equator  with 
the  appropriate  frequency  (UMS)  and  electron  energy,  the  initial 
pitch  angle  would  have  to  be  large  (which  leads  to  an  unrealisti¬ 
cally  large  initial  energy  and  energy  loss).  For  L=4.48,  Fig.  7 
shows  that  resonance  with  the  observed  energies  can  occur  in  the 
vicinity  of  6000  km. 

The  equatorial  region  is  usually  selected  as  the  location  of 
a  possible  wave-particle  interaction  region  on  the  basis  that  the 
resonance  conditions  remain  relatively  constant  for  a  greater 
distance  along  the  field  line  there  than  elsewhere,  since  B 
(which  determines  the  local  cold  electron  gyrofrequency)  and 
are  relatively  constant  there  and  a  and  6  change  more  slowly 
there.  The  L=1.66  curve  in  Fig.  7  indicates  that  for  a  small  re¬ 
gion  around  the  equator,  the  resonant  energy  decreases  slightly 
as  the  particle  goes  away  from  the  equator,  enabling  the  particle 
to  stay  in  resonance  as  it  loses  energy.  Evaluation  of  precipita¬ 
tion  events  in  the  inner  zone  (Koons  et  al  [1981])  indicated  that 
the  resonance  region  was  near-equatorial.  This  is  also  true  of 
low  energy  electron  precipitation  in  the  lower  portion  of  the 
slot  [Imihof  et  al,  1983]  (which  may  have  been  just  a  continuation 
of  the  inner-zone  precipitation,  extending  to  very  low  energy  and 
therefore  up  to  slot  altitudes) .  But,  for  energetic  electron  pre¬ 
cipitation  in  the  central  portion  of  the  slot,  E>200  keV,  L~2.3, 
resonance  conditions  can  be  met  at  the  equator  only  for  a  high 
order  resonance  or  unrealistic  conditions.  However,  examination 
cf  the  resonance  conditions  shows  that  an  electron  going  down  the 
field  line  may  continue  in  resonance  as  long  as  the  rate  at  which 
it  is  losing  energy  (which  decreases  a  and  P)  just  balances  the 
changes  in  6,  which  decreases,  and  |i.  and  both  of  which  in¬ 

crease  at  lower  position  on  the  field  line.  Since  the  potential 
well  of  the  wave  may  also  be  converting  transverse  energy  to  par¬ 
allel  energy  elastically,  the  picture  can  become  quite  compli¬ 
cated,  A  thorough  discussion  of  this  point  is  beyond  the  scope  of 
this  paper,  but  the  observation  of  a  Siple  Station  VLF  transmis¬ 
sion  being  amplified  prior  to  its  first  passage  through  the  equa¬ 
torial  region  [Inan  et  al,  1977]  can  be  taken  as  evidence  that 
off-equatorial  wave-particle  resonant  interactions  do  occur. 


m 


DISCUSSION  and  SUMMARY 


* 


Pv'- 

%  •y 


v>. 

' » 

V'v’ 

y*/’’ 

V'V' 


Drift-rate  analysis  of  electrons  observed  in  the  drift  loss 
cone,  coupled  with  the  operating  sequence  of  the  VLF  transmitter 
UMS,  explains  the  detailed  electron  observations  as  a  function  of 
energy  for  several  different  L-shells.  Additionally,  the  pitch- 
angle  distributions  which  are  observed  are  consistent  with  those 
that  would  be  produced  if  the  electrons  had  been  precipitated  in 
the  vicinity  of  UMS.  The  combination  of  pitch-angle  and  energy 
observations  are  consistent  in  every  detail  with  a  scenario  in 
which  UMS  precipitated  electrons  over  a  wide  energy  range  in  sev¬ 
eral  different  L-shells,  including  the  outer  zone. 

A  number  of  questions  remain:  a)  Why  has  precipitation  of 
outer-zone  electrons  by  VLF  transmitters  not  been  observed  previ¬ 
ously?  b)  Why  does  the  precipitation  at  a  given  energy  occur  at 
discrete  L-values  instead  of  being  a  continuum  over  the  slot  and 
outer  zone?  Is  this  related  to  the  plasmaspheric  structure?  c) 
How  important  is  the  precipitation  of  electrons  by  VLF  stations 
in  the  outer  zone  and,  more  importantly,  in  the  morphology  of  the 
m.agnetospheric  particle  population?  Is  it  the  primary  controlling 
factor? 

The  first  of  these  questions  may  have  the  simplest  answer: 
Outer  zone  precipitation  of  electrons  by  ground-based  VLF  trans¬ 
mitters  may  not  have  been  observed  previously  because  almost  no 
one  has  looked  for  it.  The  SEEP  experiment  (Imhof  et  al  [1983]) 
attem.pted  to  observe  outer  zone  precipitation  and  was  unsuccess¬ 
ful,  (The  reason  for  this  lack  of  success  might  be  because  NAA 
does  not  precipitate  outer  zone  electrons  as  does  UMS  or  it  might 
be  due  to  a  lack  of  sufficient  intensity  in  precipitated  outer 
zone  electrons  during  that  campaign  to  be  observed  by  the  SEEP 
experiment.  In  the  S3-3  data  set,  only  7%  of  the  orbits  had  outer 
zone  precipitations  which  were  thought  to  be  due  to  VLF  transmit¬ 
ters.)  In  the  outer  zone,  a  number  of  other  processes  precipitate 
electrons.  During  magnetic  storms,  intense  precipitation  at  all 
energies  is  observed  at  high  latitudes.  During  minor  distur¬ 
bances,  the  outer  edge  of  the  outer  zone,  as  observed  by  low  al¬ 
titude  satellites,  moves  in  and  out  significantly  due  to  rela¬ 
tively  low  altitude  processes  [Vampola,  1977b].  In  order  to  es¬ 
tablish  that  a  given  precipitation  event  is  due  to  a  specific  VLF 
transmitter,  a  lengthy  pitch-angle  distribution  analysis  must  be 
performed  w.hich  is  not  easily  automated  because  it  must  identify 
the  structure  as  being  in  the  drift  loss  cone,  as  having  a  dis¬ 
tribution  that  is  narrower  than  the  local  loss  cone,  and  as 


31 


was  not  de  t  e  r  rr.  1  n  e  c 


the  Sc 


lar.tic  Ar.crr.aly  (cr  a  similar  magnetic  structure  in  the  vicinity 
of  160-  El).  Most  electron  spectrometers  flown  on  satellites  do 
net  obtain  complete  pitch  angle  distributions.  For  measurem.ents 
miade  at  higher  altitudes  in  the  stable  trapping  region,  high  an¬ 
gular  definitio.n  is  required  in  order  to  utilize  the  trace-back 
technique.  Even  with  an  appropriate  data  set,  without  the  suspi¬ 
cion  that  a  given  precipitation  event  had  as  its  source  a  VLF 
transmission,  the  extensive  analysis  required  to  establish  it  as 
such  would  not  be  done.  The  present  analysis  is  the  result  of  a 
survey  of  data  looking  for  slot  precipitation  events  due  to  UMS . 
In  the  course  of  the  survey,  it  became  evident  that  outer  zone 
precipitation  in  the  few-hundred  keV  range  was  frequently  related 
to  the  presence  or  absence  of  an  energy-dependent  slot  peak.  Fur- 
therm, ore,  the  outer  zone  peak,  while  not  appearing  to  be  L-depen- 
dent,  exhibited  strange  energy  behavior  (occasionally  low  energy 
electrons  were  not  observed  and  in  some  events  the  433  keV  chan¬ 
nel  had  miuch  higher  fluxes  than  the  235  keV  channel)  . 
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The  second  question  is  more  difficult  to  answer.  The  ray¬ 
ing  "exercises"  were  instructive  and  offer  several  possible 
ers.  First,  under  som.e  conditions  of  injection  latitude,  in¬ 
ion  angle,  and  plasmaspheric  configuration,  the  waves  follow 
field  line  quite  closely  for  a  significant  distance  and  pro- 
a  large  path-length  and  interaction  time  with  coherent  con- 
ens  and  very  slow  change  in  resonant  energy  which  may  enhance 
tering  for  low  signal  levels.  Perhaps  waves  must  be  injected 
specific  latitude  and  plasmaspheric  conditions  to  produce 
ipitaticn.  Then  injection  of  waves  at  various  latitudes  would 
esponsible  fer  the  several  L-regions  of  precipitation.  A  sec- 
observation  was  that  waves  injected  at  several  latitudes  may 
pass  through  the  same  region  of  the  magnetosphere,  a 
using"  effect,  which  could  provide  constructive  interference 
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ion  occurs  on  discrete  L-shells  due  to  discrete  regions  of 
tructive  interference.  The  loci  at  which  the  several  ray 
s  intersect  are  at  high  latitude  well  away  from  the  equator, 
her  possibility  is  that  waves  aniplified  by  interacting  with 
trons  at  low  L  (e.g.,  inner  zone,  equator)  are  later  re¬ 
ted  at  high  latitude  and  travel  back  upward  toward  the  equa¬ 
te  again  interact  with  electrons  on  a  higher  L-shell  (but  not 
ssarily  near  the  equator).  These  waves  could  be  again  ampli- 
,  again  reflected,  and  again  interact  with  electrons  on  a  yet 
er  L-shell.  Assuming  sufficient  reflection  (the  product  of 
ification  and  reflectivity  being  greater  than  unity),  the  in- 
ctitn  could  be  between  upward-going  waves  and  downward-going 
ides,  as  assumed  for  the  earlier  discussion,  or  with  upward- 
el  dr.g  particles  a.nd  downward-going  waves,  in  either  hem.i- 


sphere.  Thus  the  pattern  of  precipitation  in  L  could  be  due  to 
the  specifics  of  magnetospher ic  wave  propagation  from  a  single 
injection  point  or  multiple  injection  points.  Inhomogeneities  in 
the  ionosphere  and  plasmasphere  could  play  a  crucial  role.  Fi¬ 
nally,  it  may  be  possible  that  electrons  are  indeed  precipitated 
throughout  the  magnetosphere  above  L=2  but  their  appearance  in 
the  drift  loss  cone  is  dependent  on  the  intensity  at  higher  alti¬ 
tudes.  Further  investigations  along  this  line  are  planned  but  the 
present  ray-tracing  codes  will  have  to  be  integrated  with  real 
miagnetic  field  models  and  in-situ  measurements  of  electron  den¬ 
sity  profiles  in  order  to  provide  reliable  numerical  results. 

The  last  question,  how  important  is  the  precipitation  of 
outer-zone  electrons  by  VLF  stations,  is  probably  answerable  by  a 
numlDer  of  means:  1.  Turn  off  all  VLF  stations  for  a  few  months 
and  observe  the  behavior  of  the  energetic  electrons  in  the  magne¬ 
tosphere.  (Use  only  satellites  for  navigation  and  communication 
during  that  period.)  This  method  has  the  disadvantage  that  if  the 
VLF  transmitters  are  crucial  in  determining  the  structure  of  the 
electron  belts,  the  inner  zone  could  become  extremely  hazardous 
to  satellites  through  a  build-up  of  relativistic  electrons  dif¬ 
fusing  inward  from  the  outer  zone.  2.  Do  an  inventory  of  elec¬ 
trons  in  the  drift  and  bounce  loss  cones  to  determine  where  they 
left  the  stable  trapping  region.  (Previous  studies  have  tended  to 
look  for  the  longitude  at  which  the  drifting  electrons  would  en¬ 
ter  the  atmiosphere,  since  that  is  a  much  more  easily  determined 
parameter.)  This  is  a  large  undertaking,  but  would  not  require 
the  cooperation  of  the  various  nations  operating  high-powered  VLF 
transmitters  and  would  not  be  potentially  hazardous.  3.  Make  si- 
m.ultaneous  measurements  of  high  altitude  and  drift  shell  electron 
fluxes  during  a  magnetically  quiet  period  and  attempt  to  deter¬ 
mine  the  lifetime  of  electrons  due  to  removal  by  VLF  transmission 
precipitation  and  compare  those  lifetimes  to  radial  diffusion  co¬ 
efficients  determined  from  the  same  data  set.  This  study  would 
require  a  dense  data  set,  probably  obtained  simultaneously  by 
several  satellites. 
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T^'»'  Aer'sp.rikf  i^irpor.^c  i  *'1  r 's^'u- 1  i -'fis  as  in  ar  •  h  i  r c  i -n’t*  r  tor 

:'.atl  nil  siH'urltv  pr.>ierts,  sper  i  a  1 1  ^  i  in  aivan.-eti  'T^illtary  space  'iwtems. 
•r^vijme  research  supp-irl,  the  c..rp')rat  ion  ’  s  l,ah<^  ri  t  or  v  i>perati  'ns  'n.lucts 
e  X  pe  r  1  net  t  a  1  itvl  the*>retlca!  i  nves  t  i  g  <  t  i  )ns  that  tocus  <ji^  the  app  1  *  1 1 of 

;  1  e'nc  1  M  ■'  anii  techni  H'.’^nees  to  sij<.  h  systems.  Vit<l  to  the  success  ut 

these  invest  i^at  Ions  is  the  technical  staff's  wide -ran^.  1  ni4  eK-pertisi-  and  its 
ahiLitv  tv>  Slav  ctirrent  wi  t  ti  new  developments.  This  expertise  Is  enha:i>'-‘d  Sv 
a  research  program  aimed  at  dealing  with  the  many  prohlems  associated  with 
rnpidlv  evi)lvin.t  space  systems.  fJont  r  ihut  i  their  apahlltties  t  >  t  tie 
research  eM  >rt  are  these  individual  laboratories: 

\e  rophvs  1  c  s  Lahoratory  ■  Launch  vehicle  and  reentry  fluid  mechariics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dvnamics,  env i r onmenta I  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  con t amlnat i on ,  thermal  and  structural 
control,  hij^h  temperature  thermoraechanlcs ,  kinetics  and  radiation;  cw  atvi 

pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  c>>ncr')l,  atmos'pheric  prufMgat  i  on ,  ^^ser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state -specif tc  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f  ield -ot  -view  rejecti/)n, 
applied  Laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  bittery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
env ironmencai  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  appl  icatlons,  communlcat ion  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microivave  semlcondoctor  devices,  mJcrowave/njl  1 1  imeter  wave  measurements, 
diagnostics  and  radloraetry,  mlcrowave/ml 1 1 imeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratorj^:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy -Induced 
env  i  r'jnment  s . 


Space  Sciences  l..aboratory :  Magnetospher  Ic ,  auroral  and  cosmic  ray 
pfwsics,  wave -part  1 c le  Interactions,  magnetospheric  plasma  waves;  atmospheric 
and  i'^nospherlc  physics,  density  and  composition  of  the  u;>per  atmosphere, 
remote  sensing  using  atmr>spher -c  ra<1latlon;  solar  physics,  infrared  astronumv, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  stnrms  and 
nui'lear  explosions  on  the  earth’s  atmosphere,  Ionosphere  and  magnetosphere; 
effects  »f  e  1  ,>cC  rtmagneC  i  c  and  particulate  radiations  on  spare  systems;  space 
i ns t r umen t at  ion. 
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